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Abstract 

We expand on previous work involving "vacuum-bounded" states, i.e., states 
such that every measurement performed outside a specified interior region 
gives the same result as in the vacuum. We improve our previous techniques 
by removing the need for a finite outside region in numerical calculations. We 
apply these techniques to the limit of very low energies and show that the 
entropy of a vacuum-bounded state can be much higher than that of a rigid 
box state with the same energy. For a fixed E we let L[ n be the length of a 
rigid box which gives the same entropy as a vacuum-bounded state of length 
L m . In the E — > limit we conjecture that the ratio L' m /L m grows without 
bound and support this conjecture with numerical computations. 
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I. INTRODUCTION 

In a previous paper we defined a "vacuum-bounded state" to be a generalized state 
(i.e., density matrix) for which every operator composed of fields at points outside a specified 
interior region has the same expectation value as in the vacuum. We used a vacuum-bounded 
state as a model of the state resulting from the complete evaporation of a black hole. We 
took the interior region to be a sphere of radius R, the distance that information might have 
propagated after unknown physics came into play in the evaporation process. 

In the present paper we redo our analysis of maximum-entropy vacuum-bounded states 
in a way that directly handles an infinite outside region without the need to put the entire 
system in a box. Using this technique we calculate the maximum entropies of vacuum- 
bounded states with very low energies. As in we do our calculations with a single scalar 
field in one dimension. We argue that as the temperature of the system becomes very low, 
the entropy is much larger than what one would have in a system with a rigid boundary, 
and that this difference grows without bound as the temperature goes to zero. 

As in we approximate the continuum with a lattice of harmonic oscillators. Even 
with an infinite vacuum outside our specified region, we can reduce the problem to a finite 
number of degrees of freedom. We then compute the solution numerically to support, in 
this approximation, our claims of high entropy at low energies. 

We also use these techniques to recalculate the bounds we derived in 0. We obtain the 
same bound as in that paper, but now we work directly with an infinite outside vacuum, 
which we could do only by extrapolation in |jj . 



II. PRELIMINARIES 



We work with a one-dimensional scalar field with classical Hamiltonian 

/ 7 I \ 2" 
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(2.1) 



We let L; n be the size of the inside region. We will initially put the system in a box of length 
L, so the outside region stretches from x = L UI to x = L, but we will shortly take L to oo. 

Our goal is to find a density matrix p such that if out is any operator composed of fields 
at points Xj > L in , then 



Trp0 out = (0|O O ut|0). 
From we expect p to be a Gaussian of the form 



(2.2) 



(2.3) 



where ranges over the quadratic operators <p(x)<p(y) and 7r(x)ir(y) with x,y > L m . The 
coefficients f a must then be chosen so that p satisfies Eq. (p.2[). 
In discrete form, we have [l.j 



H=-(P x -P x + x-Kx) 



(2.4) 
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where g = 1/ L\ and L\ = L/(N + 1) is the lattice spacing. 

We will take Ni a ~ L- m jL\ of the oscillators to represent the inside region, and N Q 
N — N in to represent the outside region. Equation ( |2.2|) then becomes 



Tr pXiXj = (0\xiXj\0) 
Tr pPiPj = (0\PiPj\0) 



(2.6a) 
(2.6b) 



where i and j run over the oscillators which represent the outside region. We require also 
that p satisfy an energy constraint, 

Tr pH = E . (2.7) 
We will define matrices X and P whose elements are the quadratic operators via 



= P P 
1 U 1 V 



p 



so that Eqs. ( |2.6a| ) and ( |2.6b| ) become 



Tr pX out 

,out 

Tr pP ou t,out 



(0|X out)OUt |0) 

<0|P O ut,out|0). 



(2.8a) 
(2.8b) 



(2.9a) 
(2.9b) 



We expect p to have the form 



p oc e 



-0{H+fijXiXj+gijPiPj) 



(2.10) 



III. REDUCTION 

Equation ( |2.10|) has A^ out (A^ out + 1) degrees of freedom in the / and g parameters. We 
would like to reduce these degrees of freedom to a number which does not depend on iV out , 
but only on iV in , so that we can take N out to infinity and still have a well-defined problem. 
To make this reduction we will describe our system by a new set of modes, such that the 
number of modes which can be excited in any p which satisfies Eq. ( |2.9| ) depends only on 
N- m . 
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A. A different description of the vacuum 

We will start by looking at the vacuum state. We can make a change of coordinate to 
put the Hamiltonian of Eq. fl2.4|) in diagonal form. Let Z be a matrix whose columns are 
the eigenvectors of K, K = ZQ^Z^ 1 , with the normalization 

Ztt Z T = I and Z~ lT ti^Z' 1 = K . (3.1) 

T 

Define new coordinates z via x = Zz and P x = Z^ 1 ~P Z . In these coordinates, 

H=lj2^ ] {Pi+zl)- (3-2) 

The vacuum is the ground state of this Hamiltonian. We can define raising and lowering 
operators 

a a = (z a + iP Za ) (3.3a) 

ai = ±={z a -iP Za ) (3.3b) 

# = E4 0) (aW£) • (3.3c) 

The vacuum is the state |0) annihilated by all the a a . It is straightforward to write the 
expectation values in the vacuum state, 

(Z) Q/3 = (0\z a z p \0) = ^5 a(3 (3.4a) 

(F z ) a p = (0\P Za P Zf3 \0) = ^5 aP (3.4b) 

so 

(0|X|0) = Z(0\Z\0)Z T = l -ZZ T (3.5a) 

(0|P|0) = Z~ xT 1$\$> Z \$)Z- X = l -Z- xT Z- x = ^((OIXIO))- 1 . (3.5b) 

There can be many different Hamiltonians that have the same ground state. If we 
consider 

H' = ^(P x -T'P x + x-K'x) (3.6) 

with T' and K' some coupling matrices, we can follow the above derivation to get a normal 
mode matrix Y and some frequencies Q with 

(3.7a) 
(3.7b) 
(3.7c) 
(3.7d) 

% + (3-7e) 
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= K' 
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We can define raising and lowering operators for these modes, 



1 

6/3 = 71 ( yf3 + iPy ^ 


(3.8a) 


ha = —p (ye - iPy*) 
v2 ^ ' 


(3.8b) 




(3.8c) 


= W _1 z where 






(3.9) 



Similarly, P y = Y T P X = Y T Z~ l P z = W T P Z . Consequently we can define a Bogoliubov 
transformation, 



= ^(^K + 4) + ^ 
= \ ((W- 1 + W^)^ + (W- 1 - W T ) Pa ai) . (3.10) 



For H and if' to have the same vacuum we require that the bp depend only on the a a and 
not on the a' a , which is to say that W~ x = W T , i.e. that W is a unitary matrix. With W 
unitary, YY T = ZZ T so the vacuum expectation values of Eq. ( |3.5|) have the same values 
expressed in terms of Y as they had in terms of Z. 



B. Modes that remain in the ground state 



Now consider a unitary matrix W and let Y = ZW as in the last section. Suppose we 
can find W such that Y has the following property: 

The N modes can be divided into N gs > "ground state" modes and iV free = 
N — N gs "free" modes such that for all a < N in and for all (3 > N {ree , Y af3 = 
and Y^ = 0. 

That is to say Y and Y^ 1 will have the form 

iVfr, 

( ~\ 

Y ~- 



and 



Y 



-i 



^in,free 





^out,free 


^OUtjgS 


Win 

free, in 


Nout 

v-i 

free, out 


v 


Y- 1 

± gs,out 



j}iVin 



(3.11) 



'}N S s 



(3.12) 
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If this is the case, then for (3 > N {rcc , bp = \Y^x 1 + iY y pP x ^j / y/2 depends only on 
outside operators Xi and P Xi . In any vacuum-bounded state, regardless of entropy consid- 
erations, (xiXj) and (PiPj) have their vacuum values. Consequently, for all (3 > iVf rcc , if p 
describes a vacuum-bounded state then Tr pbpbp = (0 1 fc»j 3 £» / 3 1 0} = 0. The vacuum-bounded 
constraint forces modes (3 > N {ree to be in their ground states. These modes will not con- 
tribute to the calculation of the maximum-entropy vacuum-bounded state. 

How many such modes can exist? Let Wp denote a column of W with (3 > Nf ree and let 
Z a denote a row of Z with a < N- m . Similarly let (Z~ 1 ) a denote a column of Z^ 1 . Since 
Y = ZW, Y a p = whenever Wp is orthogonal to Z a . Similarly Y^ 1 = W T Z~ l so Yg" 1 = 
whenever Wp is orthogonal to (Z~ l ) a . Since W is unitary, the Wp must also be orthogonal 
to each other. Thus there are N gs columns of W which have to be orthogonal to N- m rows of 
Z, to N in columns of Z~ l , and to each other. Since there are N components in a column of 
W it can be orthogonal in general to at most N — 1 other vectors. Thus iVg S is limited by 
N — 1 — 2N in + N gs - 1, or N gs = N - 2N in = N ont - N in . Thus whenever iV out > iV in there 
will be iV gs = iV ou t — N- m modes that are forced to remain in the ground state. 

These conditions determine the columns Wp with (3 > 2N in up to a unitary transforma- 
tion on these columns alone, and likewise the remaining columns are determined up to a 
unitary matrix which combines them. 

C. Density matrix and entropy 

Let Tt be the Hilbert space of states of our system. We can write TC = 7i <g> Q where 7i 
is the Hilbert space of states of the "free" modes and Q is the Hilbert space of states of the 
"ground state" modes. Let |0) gs G Q denote the ground state of system Q. For any operator 
A we can define an operator A that acts on H via (<S|A|/?) = (a <g> Q gs \A\/3 ® gs ) for all 

a, pen. 

Let p describe a vacuum-bounded state and let p gs = Tr frcc p. Then Tr p gs bpbp = 0. This 
defines the vacuum state, so p gs = |0) gs (0| gs and thus 

P = P®(|0) gs (0| gs ) . (3.13) 

If we write p in diagonal form, 

p = J2P«\a)(a\, (3.14) 

a 

then 

p = ^P a |5® gs ) (a® gs | . (3.15) 

a 

The entropy is 

S = -Trplnp = -^s(Pa) = -Trplnp. (3.16) 

a 

Since any vacuum-bounded p has this form, a variation of p that preserves the constraints 
must also have this form, 
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6p = 6p® (|0) gs (0| gs ) . (3.17) 
If S = -Tr (p + tip) ln(p + tip) = -Tr (p + tip) \n(p + tip), then 

^- = -Ti6plnp. (3.18) 

If p maximizes S subject to the constraints, then we must have Tr^plnp for any variation 
5p that preserves the constraints, i.e. for which Tr 5p = and 

Tr 5pH = Tr bpx i x j = Tr bpP { Pj = . (3.19) 

Now for any A, Ty 5pA = Tr (Sp ® (|0 gs )(0 gs |))A When we take the trace we only need to 
sum over states of the form \at ® gs ). Thus 

(A) = TidpA = ^(a\8p\(3)0® O gs \A\a® gs ) = Tr bp A = (A) . (3.20) 

a/3 

Thus we are looking for p that maximizes S = — Trplnp subject to the constraints 
Tr SpH = Tr SpXjXj = Tr 5pPiPj = 0, where i and j range over the outside oscillators. 
From Eq. ( |2.10|) we expect p to have the form 



p oc e V / . (3.21, 

As in |l]] we can write this 

p oc e- fia ' , (3.22) 
where H' is a fictitious Hamiltonian for these oscillators, 

H' — H + fijXiXj + gijPjPj ■ (3.23) 



D. New coordinates 



We would now like to introduce new coordinates w as follows: The first iVi n w coordinates 
will be the inside oscillator coordinates, w in = x in . The last iV gs coordinates are the ground 
state normal modes, w gs = y gs . The remaining N in coordinates can be any coordinates that 
are independent of those specified so far; we will make a particular choice later. 

To do this, we proceed as follows: Using Eqs. (|3.11|) and (|3.12|) we can write out 



y— 1 y: _i_ -y— 1 y 

1 free, in- 1 m,free T 1 frccout- 1 out, free 



V — i y 

1 gSjOut- 1 out, free 



y— 1 y 

1 free,out 1 out,gs 



y— 1 y 

- 1 gs^ut- 1 out,gs 



(3.24) 



In particular, ^ S)0 ut^out,gs 
matrices i? and with 



/. We would like to extend Yg Si out an d y"out,gs into square 
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N in N gs 

R= ( D\Y out , gs )} Nout (3.25) 



and 

N out 



JT> = (") \ N - . (3.26) 

\ 1 gS,OUt / J Jv m 



This means that we must find D and D' such that 



Y-^ t D = (3.27a) 
DT out , gs = (3.27b) 
D'D = I. (3.27c) 



There are many possible choices of D and D' that satisfy Eqs. ( [3.27D . Here we proceed as 
follows: Let 



Z m 



^out 



(3.28a) 



Z- X = (Z£\Z^) (3.28b) 
so that Z m Zr^ = 7, Z out Z~^ t = I and Z^Z- m + Z~±\Z ont = I. Now let 

D = l -Z out Zl = (0|X|0) out)in (3.29a) 

D' = \zrfz-^ t = (0|P|0) in>out , (3.29b) 

let A and B be N- m x N- m matrices with 

AB = (D'D)- 1 (3.30) 

and let 

D = DA (3.31a) 

D' = BD'. (3.31b) 

We also divide 

W = (W iree \W gs ) . (3.32) 
Since Y = ZW and Y~ l = W T Z~\ we have 



Z in W gs = (3.33a) 
7 T Zr l 

gs in 



WlZ^ 1 = . (3.33b) 
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Thus 



D'Y mt>ss oc BZrf Z^ t Z ont W ? « = BZr 1 



in 



J out Z/ ouf>'gs 

by Eqs. Q3.33| ) and their transposes. Similarly 



BZr' Zr^Wv = 



Y-'D oc WlZ~lZ mt ZlA = WlzlA - WlZr}Z- m Z?A = . 



(3.34) 



(3.35) 



From Eq. (|3.30|) we find D'D = I. Thus the matrices D and D' satisfy Eqs. ( [3.271) . We still 
have the freedom of choosing the matrix A arbitrarily. 
Now let 



Q 





N out 










R 



(3.36) 



1 

and define w by x = Qw so P x = Q~ x ~P W . Then w in = x in and w gs = y gs as desired. 



E. Reduced operators 

We would like to recast our problem in terms of w free , the the first iVf ree w coordinates. 
First we look at the operators x^x u = X^ and P^P V = P^. If we write 



Q = (<2free|Qgs) (3.37a) 

Q- 1 = (^f) , (3.37b) 
we have x = Qw = Q ivce w fvee + Q gs w gs and P = Q~ lT ~P w = <3^J( p ™)free + Qg S 1 (Pw) gs , so 



X = QfreeWfree^eeQfree + Qg S W gSjgs Q gs (3.38a) 

? gs (Pw)gs,gsQ gs 



1^ 1 1^ 1 

^ = Qfrcc (P«j)free,freeQf ree + Q RS (^w) gs,gsQ gs ■ (3.38b) 



Now we form the reduced operators x~^x v and P^P U . Since w gs = y gs and the y gs modes are 
in the ground state by definition, W gSjgs = (Pu,) gsgs = (1/2)/, and so 

X = QfrceW freeifree Q£ ee + ^Q g sQ gs (3.39a) 

~ 1 1^ 1 

^ = Qfree O^w) free,freeQf ree + gs ' (3.39b) 

In each case the second term is just a constant. Now 

H=~ {P,P, + K^x^Xv) = ^Tr (P + KX) (3.40) 
where the trace is over the oscillator indices. Thus 
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where 



H — -Tr ( Q fr c C \ Li; »• 'Irci JVcc 



Q + A"Q frcc W frcc>frcc (5fr ec ^ + const 
^Tr (r(P w ) freei f r ee + ^W frcc , frcc ) + const 



T QhccQiTCC ^ g 



N in 


Win 










D'D lT 



(3.41) 



(3.42a) 



A — Qf ree KQf ree 
The constant term in Eq. Q3.41 ) is 



- 11 -in,in 




^ - I1 out,m 


D T K outout D 



^Tr ( KQ gs Q T gs + g gs lJ g- 1 



— Tr I V T ft" V _|_ V -1 V -1 



out,gs 



gs,out gs,out 



(3.42b) 



(3.43) 



It depends on the ground state modes only and is just part of the zero-point energy. It will 
be the same in the vacuum and in a vacuum-bounded state. Thus if instead of Eq. ( |3.41| ) 
we use 



(3.44) 



H = -Tr (T(P tt ) fr ee,iree + AWf re e,free 



we are just shifting H by a constant term and thus changing the zero-point energy. 

Now we would like to make this reduced system look as much as possible like the system 
we started with. Let Bi be some matrix such that B\B\ = (D'D' T )~ l , let B 2 be a unitary 
matrix to be determined, and let B = B 2 B\. Then D'D' T = I, so T = I. This gives 
A = A X A 2 where A 1 = (D'D^B^ 1 and A 2 = Bj. Let D 1 = DA 1 and let 



m,in 



T) T K 



Aj nou t-Dl 



Now 



in, out 



-9\ 



(3.45) 



(3.46) 



so K in ont Di is nonzero only in the last row. Thus the last N in + 1 rows and columns of K\ 
look like 



1 A^in 



2<? 


? 


? 






(3.47) 



A matrix of size (N in + 1) x (N in + 1) can be put in tridiagonal form by a unitary transform 
of the form 
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(I 






V 





(3.48) 



which can be constructed, for example, using the Householder process. We will use this to 
choose A\ = U, so that K will be tridiagonal. For each off-diagonal element in the resulting 
tridiagonal matrix there is a choice of sign, and we will choose them all to be negative. Thus 
in the w coordinates, each oscillator has unit mass and is coupled only to its neighbors. 
This completely specifies the matrices A and B and thus D and D'. Note that D and D' do 
not depend on the undetermined parts of W. 

To make H' in Eq. (|3.23| ) we can add to H a kinetic and potential term involving outside 
oscillators only. The potential term (disregarding a constant) is 



fijXiXj — Tr /X c 



Tr/Q 



out, free free, freeQout free 



(3.49) 



Since 



Win N in 



<2out,frec = ( | -D )}iV 



(3.50) 



this term is equivalent to adding an arbitrary term to just the lower right part of K. Similarly 
the kinetic term is 



g ij P i P j = TrgF, 



out, out 



1 1 

Tr #Qfree,out(^w)free,freeQfree. 



out 



(3.51) 



Since 



Nor, 



Qftet 



\ }n ou 



free, out 



D' J }N in 



(3.52) 



this term corresponds to adding an arbitrary term to just the lower right part of T. 
That is to say we can write 



H' = -Tr (T'(Pj frceifree + iTW f ree,free) 



(3.53) 



with 



r 

K' 














-^in,in 


-Kin,out-D 


-D^i^out.in 


? 



(3.54a) 
(3.54b) 
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F. Reduced constraints 

Now we rewrite our constraints, Eqs. ( |2.6| ), in terms of the w coordinates. We will keep 
only the parts of the constraint equations that are not automatically satisfied because of the 
ground-state modes. For the expectation value constraints, from Eqs. ( |3.2(J| ) and ( |3.39| ) we 
have 

(X) 

out, out (^0 out, out Q out, free (^0free,freeQ O ut, free ^ ^ out iS s ^out,gs 

= D(W) mid , mid D T + const (3.55a) 

2 

= D' T (F w ) midimid D' + const (3.55b) 

where w mid means the outside elements of w free , i.e. Wjv in +i • • • w 2N in - 

These expectation value matrices must be the same in the vacuum as in the vacuum- 
bounded state. We can accomplish this by requiring that (P«,) m id,mid and (W) m id, m id are 
the same as in the vacuum. Our problem is now equivalent to one with only N- m outside 
oscillators, and thus only iVi n (iVi n + 1) expectation value constraints, regardless of the value 
ofAUf 

For the energy constraint, Eq. ( |2.7| ), we are concerned only with the renormalized energy 
Tr pH — (0\H |0). Thus the constant term in Eq. (|3.41|) does not contribute, and we can use 



- rp XT 

(P) 

out, out (^)out.out ^ free, out (P?i>) free, free^free, out r\ ^gs,out^5gs,out 



H from Eq. ( |3.44| ). Once again there is no dependence on N Q 



G. Derivation based on inside functions 

We would now like to redo the proceeding calculation in a way which does not depend 
on the number of outside oscillators. Then we can remove the box from our system by 
taking L — > oo and N — > oo. It appears that we have used the matrices D and D' which 
have an index that runs from 1 to iV out . However, we have used them only in particular 
combinations. The quantities which we need in our calculation are 

1- -^in, out -D 

2. D'D 

3. D'D' T 

4. D T K onttOUt D. 

Each of these quantities is an iVi n x iVi n matrix, so it is reasonable to imagine that they do 
not depend on iV out in the N — > oo limit. 

We proceed as follows: From Eqs. ( |3.29| ) we have 

D = (X) ouMn (3.56a) 
D' = (P) in , out • (3.56b) 

Any given element of (X) and (P) has a smooth limit when iV out is taken to infinity. It is just 
a particular expectation value of a half-line of coupled oscillators, which is a well-defined 
problem. We can express the above items in terms of such elements as follows: 
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1- -K'in.out-D depends only on the first row of D which is (X)jv in +i,i n so it is already well- 
defined in the limit. 



2. From Eq. ( |3.5b| ) we have 



>)<X> = il (3.57) 



so D'D = (P)in,out(X)out,in = \l ~ (P)in,m(X)i n ,i n , which does not depend on iV out . 



iT 



3. We expand (P)(P) = \Z~ l Z^Z' 1 Z~\ We insert / = ZVt Q Z T here to get 

(P)(P) = l -Z- lT n,Z- l = - A K. (3.58) 

Then D'D' T = (P) in , ut(P)out I i n = \K inM - (P) in , in (P) in , in which does not depend on 

iVout- 



4. Using Eqs. ( |3.57|) and ( |3.58| ) we can write 



so 



(X)K = K(X) = (P) and (X)if(X) = ij (3.59) 



+ (X) 

in , in -Kin , out ( X ) ou t , i 



4 



~^(^)in,out^out,in(X)i ni i n -|- {^^} in,in^"in,in (-^^} in,in 
D T K outout D + (X) i ni i n ( (P) in 5 i n — ^in,in (^) in,in) 

,in ~~ (X/in^n-Kin^n) (X)j nj i n + (X)i n> i n i^i n i n (X)i ni i n (3.60) 



and so 



D T K on t,ontD — —I — (X)i n) j n (P)i n; i n — (P)in,in(X)i n) i n + (X) i n ,in-f^in,in (X) i n i n (3.61) 



which does not depend on N, 



out • 



Thus we can now take iV — > oo with N[ n fixed and have a well-defined problem in terms 
of K with a finite number of free parameters. 

H. Calculation of the reduced vacuum 

We are trying to compute K in Eq. ( |3.42b| ) in the N ont — > oo limit. We will keep iV in 
and the oscillator spacing Li = L/(N + 1) fixed. With the normalization in Eq. (|3.1|) the 
vacuum normal mode matrix is given by 
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J /J,U 



with 



and 



Thus 



I 2 sin k u fji 
N + l 7/UJZ 



7TU 

N+l 



2(N + 1) k v 2 k v 



(X) 



Li 



N 



nv N + 



t5 



sin kafism k a v 



a=l 



smkafismkais 
2(N + 1) ^ Bin(A; a /2) 



Now we use 



to write 



cos(6> — (h) — cos(6> + 



2 sin 9 sin i 



i , Li cos fc Q (/i - i/) - cos /c Q (/i + v) 

x hv ~ 4(iV + l) 



a=1 sin(/c Q /2 
Using Eq. ( |3.66 ) again, for any number a we can write 

cos k a (a — 1) — cos/c a a 



and thus 



sin(fc a /2) 

cos k a {n — v) — cos k a (fi + v) 
sm{k a /2) 



2 sin k a (a- 1/2) 



2 51 sinife a (a- 1/2). 

a=fi— v+1 



(3.62) 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



(3.67) 



(3.68) 



(3.69) 



If we put Eq. (|3.69| ) into Eq. ( |3.67| ) and bring the sum over a inside the sum over a we get 
a sum that we can do, 



. 7ra(o-l/2) tt(2o-1) 

a=l 



4(iV+l) 



(3.70) 



Now we sum this over a. Since \i — v and \x + v have the same parity, the (— l) a term does 
not contribute and we get 



(X) 



/IV 



L x ^ 7r(2a - I) 
2s cot 



4(iV+l) 



a=fj,— v+1 



4(iV+l) • 



(3.71) 



14 



The sum over N has been eliminated. In the N — > oo limit, the argument of cot goes to 
zero and so we can use cotx = 1/x + 0(x -3 ) to get 



7T 



The sum can be done using special functions: 



(X) 



ft]' 



Li 
2tt 



j\) (u + v + - ) - -0 - v + ^ 



where ip is the digamma function, ip(x) = T'(x)/r(x). 

To compute (P) we write the inverse of the normal mode matrix, 



and 



(P) 



[IV 



7- 1 

jXV 



1 2fafr . , 

— — - sin k u v 
N + l M 



N 



N k, 



N 



- E u a sin k a \i sin k a v 

-L 1 



a=l 



— E sin — - shi sin k a v . 
1 ) — 1 2 



In this case the sum can be done directly. Using Eq. ( |3.66| ) we can write 



where 



/V 



k n 



Li(N + l) 
1 



E sin — cos k a X 



4L X (JV + 1) 
Again we take iV A to get 

1 



tt(2A + 1) tt(2A-1) „, 1NA 

cot -Att — - cot — - 2(-l) A 



4(iV+ 1) 



4(N + l) 



vrLx V2A + 1 2A- 1 



and 



(3.72) 



(3.73) 



(3.74) 



(3.75) 



(3.76) 



(3.77) 



(3.78) 



(P) 



itL x U(/i + v) 2 - 1 4(/i - z/) 2 - 1 



(3.79) 



Equations ( |3.73| ) and Q3-79|) give (X) and (P) in the N out — > 00 limit. Using these values in 
the procedure of sections p fD| through [III G| we can compute the matrix K numerically for 
a system with inside length L in but no outside box. 
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1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

x 



FIG. 1.: The ratio of the coupling coefficients in the reduced problem 
to what they would be in a regular problem of N{ ree oscillators, plotted 
against x = Lift, for L in = 1. 

IV. NUMERICAL COMPUTATION OF THE REDUCED VACUUM 

We have computed numerically the reduced coupling matrix K, by following the proce- 
dure of sections |III DHHTH| . This is a fairly straightforward problem in numerical analysis; 
the number of steps grows as N? n . However, in order to produce accurate results for N in > 6 
it is necessary to use very high-precision floating-point numbers. The necessary number of 
bits of mantissa in the representation appears to be about 10A^ in . The code was written in 
Lisp and executed on DEC™ Alpha™ workstations. 

The resulting matrix K is a tridiagonal matrix which gives a set of self-couplings and 
nearest-neighbor couplings for the fictitious oscillators Wf ree . We can express these couplings 
as multiples of the couplings for a regular chain of Nf ree oscillators with spacing L\. Thus 
we write the self-coupling as 

K m = -qU (&) (4-1) 
and the nearest-neighbor coupling as 

K^ +1 = 2gf^ 1/2 (^) . (4.2) 

These coupling coefficients converge rapidly to a universal form f(x) where Li/i — > x in 
the continuum limit. Some results are shown in Fig. 1. We can see that f(x) ~ 1 until 
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0.005 



0.004 



*. 0.003 



0.002 



0.001 




FIG. 2.: The coupling coefficient ratio for N- m = 175 in the region x > 1.8 
and the fit f(x) = 3.2(2 - x) A . 



x ~ 1.4 at which point it begins to fall and asymptotically approaches as x — > 2. For 
values of x near 2, /(x) is well fit by 

/(a;) = a(2-a;) 4 (4.3) 

with a ~ 3.2, as shown in Fig. 2. 

Some typical normal modes of the reduced vacuum are shown in Fig. 3. They are sine 
waves while f(x) ~ 1 and then begin to oscillate faster and faster as f(x) shrinks. At 
first the amplitude of the oscillations grows but for larger x it shrinks rapidly to zero. The 
wavenumbers in the inside region (and thus the frequencies) are smaller than we would find 
for a rigid box because most of the oscillations are in the part of the outside region where 
f(x) <C 1. In fact, as iV in — > oo we would expect the low-lying frequencies to go to zero, for 
the following reason. 

We can find the frequencies by computing the normal modes of a Hamiltonian 

H = [ L dxdy (T(x - y)n{x)n(y) + K(x - y)<t>{x)<t>{y)) , (4.4) 
Jo 

which requires solving the eigenvector equation 

[ L dydzT(x - y)K(y - z)g(z) = Xg(x) (4.5) 
Jo 

with the boundary conditions 
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0.00 



Mode 5 Freq 0.3367 
Mode25Freq 7.7131 




0.50 



1.00 



1.50 



2.00 



FIG. 3.: The 5 th and 25 th normal modes in the reduced vacuum, computed 
with iVjr, = 175. 



9(0) = 
g(L) = 0. 



(4.6a) 
(4.6b) 



Since Eq. (|4.5|) is a second-order differential equation we expect two degrees of freedom in 
the solution. However, one degree of freedom is manifestly the overall scale, which does not 
affect the boundary conditions. Since there are two boundary conditions but only one free 
parameter, we can expect to find solutions only for particular values of A. For example, for 



the usual scalar field Hamiltonian the general solution to Eq. ( [15) would be 



g(x) = csin(vA:r + 5) . 

To satisfy Eqs. ( |4.6| ) we need to choose 6 = and \/A = nir/L for some integer n. 
However, if we use 



(4.7) 



1 r2L iB 

H = d '<*> 



dx . 



(4.8) 



with f(x) — > a(2 — a;) 4 , as suggested by Fig. 2, we will get a continuum of frequencies. 
The problem is that since f(x) — > as x — > 2 the boundary condition there does not 
really constrain g(x). There can be arbitrary changes in g(x) near x = 2 and so Eq. 
( [4. 6 b| ) can always be satisfied. Since there is only one effective boundary condition and 
one effective degree of freedom, we expect to be able to find a solution for any A. Thus 
in the continuum limit there are modes with arbitrarily low frequencies. This is not an 
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FIG. 4.: The lowest frequency of the reduced vacuum and the fit 1.81x 



x 2 . 



-l 



unreasonable conclusion, since although the range of x is finite, we are using it to represent 
the infinite half-line. In the infinite vacuum there is no right-hand boundary condition, and 
there are modes of every frequency. 

This conclusion is confirmed by numerical results. In Fig. 4 we plot the lowest normal- 
mode frequency versus N- m . As shown in the figure, the frequencies are well fit by a curve 



ax 



— x 



(4.9) 



with a ~ 1.81. If this form is correct, in the N- m — > oo limit the lowest frequency goes to 
zero. 

If we go to a vacuum-bounded state we will introduce some finite temperature. We then 
expect that the non-zero temperature will increase the frequencies in such a way that there 
are only a finite number of low-lying modes and thus a finite entropy. However, in the 
limit where T — > we do expect an entropy-to-energy relation equivalent to a system with 
infinitesimal frequencies. We now consider this situation. 



V. THE LOW-ENERGY LIMIT 

In the low-energy regime we can make a first-order expansion around the vacuum. To do 
this we note that the only dependence on (3 in our equations (see |[|) is through expectation 
values that depend on coth(/5a; Q /2). For large f3 we can approximate 

coth^ « l+2e-*" . (5.1) 
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The change in coth(/3a; Q /2) is the largest for the smallest frequency, which we will call u\. 
We will ignore e~^ a for larger uj a by comparison with e"^ 1 . Thus we take 

5cot h^l = 2e-^ ac =2 e (5.2a) 

5coth^ = for«>l. (5.2b) 

Then we write 

-^mid.mid = ^~mid,mid + ^^~mid,mid (5.3a) 

-^■mid,mid = -^mid,mid + ^-^mid,mid (5.3b) 

where <5-Kmid,mid and <5T mi d,mid are 0(e). These changes give rise to 0(e) changes in U and 
the cu a , which in turn give rise to 0(e) changes in (w m w n ) and (P Wm P Wn ). 
Since overall {w m w n ) and (P Wm P Wn ) cannot change we must have 

= 5(w m w n ) = £ (- 5 ^U™UT: + ^SU ma U™ + ^U ma 5U na ) 

' UZK?e (5.4a) 



, .vac 
0J 1 



(X, . , , .vac vac 

+< c C/ vac ; 1 1 C/ vac - 1 e. (5.4b) 

We thus have iV in (A r in + 1) linear equations for N- in (N- ni + 1) unknown values of <5T mi d im id 
and 5i^mid,mid, which are readily solved. Since the inhomogeneous part of these equations is 
0(e), all the results must be 0(e) as well. In particular, the 5u a are 0(e). Now if T is very 
small as compared to all the cu a , then e will be small as compared to all the parameters of 
the problem, and so the first-order approximation will be good. For any fixed number of 
oscillators N in there will be some minimum frequency u>i, and if we take /3 C 1/wi we will 
always be in this regime. 

Now the entropy S depends only on (3 and the uj a . Since the modes are uncoupled, 

S = Y,Si(Pu a ) (5.5) 

a 

with 

Sl (Pu) = -Hl-e- pu ) + - F ^ Zi . (5.6) 

Since (3iv a ^> 1 all the terms are very small, and the cu 1 term dominates, 

S « S 1 (u 1 ) = (1 + Pwje-f*" 1 + 0(e- 2 ^) . (5.7) 

Since e drops exponentially with increasing /3, we expect that for f3 large enough, /36u>i <C 1 
so that 

S=(l + (3u™)e- puJ ™ + 0(e 2 ) . (5.8) 
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The value of S given in Eq. ( |5.8D is the one we would get from a rigid box with length 

4 n =7r/o;r c . (5.9) 

To approximate the energy, we proceed along the lines of section IX of JI[] . The direct 
calculation is made difficult by the fact that, while H' differs from H only by 0(e), we 
must subtract from both Hamiltonians a large ground-state energy. Instead we work by 
integrating on T. From F = E — TS and dF = —SdT we find 

E(T) = TS(T) - [ T S{T')dT . (5.10) 
Jo 



Now L' m depends only on u\ ac , which depends on N in but not on (3. If Eq. ( |5.8|) is valid 
for a particular N in at (3 = 1/T is it valid for T' < T and f3' = 1/T' > (3. Thus both 
S{T) and S(T') in Eq. (|5.10|) are just the entropy of a rigid box of length L\ n . Thus the 
entropy-to-energy relationship is just S(E) = S rh (L' in ; E), the entropy as a function of the 
given energy in a rigid box with length L\ n . 

For such a rigid box at very low energy we find 

E = ^—~u?e-^ (5.11) 




S = l + ln^r -s (5.12) 



and thus 



where uj\ = ir/L' in = u\ ac is the frequency of the lowest mode. 

Now for any given N in we get some uo\ ac . As discussed in section |V|, the larger N in we 
choose, the smaller u\ ac we will have. For iV in fixed we can choose (3 ^> l/u; vac and proceed 
as above to get a large value of L\ n . However, we are really interested in the continuum limit 
at fixed temperature. If we increase iV in with (3 fixed we will find that U\ (and eventually 
an arbitrary number of the uj a ) will become smaller than 1/(3. When this happens, the 
approximations of Eqs. ( |5.2| - |5.4| ) will no longer be good. 

However, we do not expect the entropy to decrease drastically in this limit. To make the 
entropy small would require making all the frequencies large. If the frequencies were large, 
the approximations we have used would again become valid. Then we could argue as before 
that the entropy should be large. It would be hard to have a consistent picture. 

Now consider the limit as T — > 0. For each T we start with some initial number of 
oscillators . We choose not too large, such that u\ ac ^> T. With this value of iV in , 
we find Z/ in °) ~ 7r/u;[ ac . We then let A^" in — > oo and we conjecture that the entropy does 
not change much, and thus in the continuum limit S(E) ~ S lh (L\^ ; E) . As we decrease 
T we can decrease the initial o;[ ac and so increase without bound. Thus we make the 
following conjecture: 

For a given energy E, let L\ n (E) be the length of a rigid box such that 
the vacuum-bounded state with energy E and length L\ n has entropy S(E) = 
S Th (L' hl (E);E). Then 

l im M^l = 00 . (5.13) 

E-+0 L in 
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FIG. 5.: The length L\ n of a rigid box that gives the same S(E) as a 
vacuum-bounded state at temperature T = 1/(3 and L m = 1.0. 

To support this conjecture numerically we turn to direct calculation of energy and entropy 
values for vacuum-bounded states at low temperature. For various fixed values of /3 — 1/T 
and for various numbers of oscillators we compute S and E and from them the equivalent 
length L; n . Some results are plotted in Fig. 5. While iV in is still small enough for the 
approximations of Eqs. (pT2|- |5.4|) to be valid, L\ n grows with N in . Once iV in has left this 
regime, it appears that L' in levels off. It is at least reasonable to believe that there is no 
further change in L\ n as N- m — ► oo. In Fig. 6 we plot the eventual level of L\ n versus (3. It 
appears that the limiting value of L' m grows nearly linearly with j3, and thus L\ n — ► oo as 
E — > as conjectured. 



VI. HIGH ENERGIES 

In this section we use the results of section JTT| to redo the computation from |TJ of the 
bound on the entropy of a high-energy vacuum-bounded state. The difference here is that 
L out and thus iV out are explicitly infinite, whereas in [l[| we were limited to iV out ~ N in . 

First we solve numerically for the maximum-entropy vacuum-bounded state. We first 



compute the reduced vacuum as in section IV and then we find the vacuum-bounded state, 



exactly as in section VII of with the vacuum couplings given by K. The results are 
shown in Fig. 7. As in [||] we find that the mode number spacing is very nearly uniform. 
Using this as an ansatz, we can repeat the calculation of section VIII of PJ. We will 
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FIG. 7.: The numerically computed wavenumbers compared with the 
best-fit line through the origin for L m = 1.0, f3 = 0.5, N in = 12. 
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work in the w coordinates, and set 



« n+ i> = « n+ i) vac (6.1) 



A. The vacuum 



In section D fU_ we computed the values of (X) . To convert to w coordinates we proceed 
as follows: From section |111D| , x = Qw and so 



^-out 



Rw 



out 



Y m 



'mid "I - -^out,gsWg S . 

In the numerical work we found that with our choice for D we got 



in, mid 







Since 



in, mid 



-^in,out D 



D 



it follows that the first row of D is (1, ... 0). Thus from Eq. ( |6.2|) we get 

X Nin+ i = W N . m+1 + (K ut,gs)l • W gs 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



where (l^,ut,gs)i denotes the first row of Y onttgs . Now 
Since the last term is non-negative, we have 

« n+ l) VaC < (^N in+ lV 



gs,gs 



gs,out^g S ,out J n ■ 



(1/2)1 and W gs , frce = 0, so 

(6.6) 



(6.7) 



This gives us an only upper bound on (w%. +1 ) vac , but that is sufficient to produce an upper 
bound on S(E), which is what we are looking for. 
From Eq. (|3.73|) we have 



(x 2 N . n+1 ) vac - (X) jv in +i,Ar in +i 



Li 
2tt 



* (2iV,„ + §) - * (i 



(6.8) 



We would like to evaluate this expression in the iV in — > oo limit with L in fixed. There is a 
prefactor of L±, which goes to zero in this limit, but that is just an artifact of the conventions 
we have used for the discrete problem, and will appear in the finite-energy vacuum-bounded 
states as well. For large x, 



ip(x) ~ lnx + 0(1/ x) 



(6.9) 



so without the prefactor there is a logarithmic divergence. We are interested in the lniVu 
term, and in the constant term, but we will ignore any terms of order 1/N in or lower. 
We use Eq. (|6.9[) and ip(l/2) = —7 — 2 In 2, where 7 is Euler's constant, to get 
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/~2 \vac _ ^1 

^+1; - 27r 



ln2iV in + 7 + 21n2 + ( — 

V in 



and so 



2tt 



ln8A^ in + 7 + 



(6.10) 



(6.11) 



B. Evenly spaced wavenumbers 

The calculation of section VIII B of U does not have any depedence on N, except in 
the factor L/N which here is L\. Thus we can use the same procedure with w in place of 
x, but the same mode function, to find (see 0, Eqs. (8.33) and (8.37)) 

. , . ^ L x ( 8iV norm sin vr A , A , l-e~ 47rT ' A \ 

Setting (vjjf + i) = (wj^ + i) vac from Eq. (|6.11|) and using the same approximations and 
in |L| we find 

IL <L in + -^- In L in T + O (I) (6.13) 



2nT \T 

just as in (l|]. This confirms our conclusion that the difference in entropy is bounded by 

5S(E)<l\nL in T^]:\nS*\ (6.14) 
6 6 

and the conclusions of |IJ remain unchanged. 



VII. DISCUSSION 

We have improved the calculations of []IJ by handling an infinite number of outside 
oscillators with only a finite number of degrees of freedom. This is possible because the 
total number of modes that can be excited depends only on the number of inside oscillators. 
Our technique is in some ways similar to representing an infinite half-line by a finite number 
of oscillators with larger and larges spaces between them. 

Using this technique we look at very low energy states which look like the vacuum except 
in a particular region. We find that, because of very low frequencies in the vacuum state, 
much more entropy can be stored in a small region with the vacuum-bounded condition 
than could be stored in such a region with a rigid boundary. To quantify the difference 
we let Li n be the size of a rigid-bounded system with the same entropy and energy as a 
vacuum-bounded system of length L in . We argue that one should expect the ratio L[J L m 
to grow without bound as the energy decreases. Numerical calculations lend support to this 
argument. 

We also redo the high-energy calculation from ]I| to reach the same conclusion with fewer 
approximations. 
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